Phase Synchronization of non-Abelian Oscillators on Small- World Networks 
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In this paper, by extending the concept of Kuramoto oscillator to the left-invariant flow on 
general Lie group, we investigate the generalized phase synchronization on networks. The analyses 
and simulations of some typical dynamical systems on Watts-Strogatz networks are given, including 
the n-dimensional torus, the identity component of 3-dimensional general linear group, the special 
unitary group, and the special orthogonal group. In all cases, the greater disorder of networks 
will predict better synchronizability, and the small-world effect ensures the global synchronization 
for sufficiently large coupling strength. The collective synchronized behaviors of many dynamical 
systems, such as the integrable systems, the two-state quantum systems and the top systems, can 
be described by the present phase synchronization frame. In addition, it is intuitive that the low- 
dimensional systems are more easily to synchronize, however, to our surprise, we found that the 
high-dimensional systems display obviously synchronized behaviors in regular networks, while these 
phenomena can not be observed in low-dimensional systems. 
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I. INTRODUCTION 



Synchronization is observed in many natural, social, 
physical and biological systems, and has found applica- 
tions in a variety of fields 0. The large number of net- 
works of coupled dynamical systems that exhibit synchro- 
nized states are subjects of great interest. As a theoreti- 
cal paradigm, the Kuramoto model is usually used to in- 
vestigate the phase synchronization among nonidentical 
oscillators 00- [1 However, some real physical systems 
that display synchronized phenomenon (e.g. the coupled 
double planar pendulums, the two-state quantum sys- 
tem, etc.) can not be properly described by the simply 
Kuramoto model. Actually, the Kuramoto oscillator is 
equal to the left-invariant flow on the simplest nontrivial 
Lie group [7(1) = {e^|C G R} (see the Refs. 0, || for 
the fundamental knowledge of Lie group). Accordingly, 
in this paper, by extending the concept of Kuramoto os- 
cillator to left-invariant flow on general Lie group, we 
investigate the generalized phase synchronization on net- 
works. In particular, the left-invariant flows correspond- 
ing to non-Abelian Lie groups are named non-Abelian 
oscillators, which can be applied on some non-Abelian 
dynamical systems like quantum systems. 

This paper is organized as follow: In section 2, the 
concept of synchronization of non-Abelian oscillators 
is introduced. In section 3, the analyses and sim- 
ulations on some typical examples are given, includ- 
ing the n-dimensional torus T n , the identity compo- 
nent of 3-dimensional general linear group GL(3, R) (i.e. 
GLo(3, R)), the special unitary group SU(2) and the spe- 
cial orthogonal group 50(3). Finally, we sum up this 
paper and discuss the relevance of phase synchronization 
of non-Abelian oscillators to the real world in section 4. 
Some mathematic remarks are given in the Appendix. 



II. SYNCHRONIZATION OF NON-ABELIAN 
OSCILLATORS 

First of all, we review the simple oscillator: 

e i(ut+e) =Ce iwt^ c=e i6 . (1) 

It is obvious that x{t) = tut + is the solution of the 
equation 



^= W , x(0) = 6, 
which relates to the dynamical system 



$(f,c) 



t 



(2) 



(3) 



on the circle S 1 = 17(1) = {e< |C G R}. 

In terms of Lie groups, S 1 = U(l) is an abelian Lie 
group with its Lie algebra iR = {iu\w G R}, and the 
map 



exp : iR — ► S , exp(iw) = e l 



(4) 



is the exponential map of S . The system (3) is also 
called the left-invariant flow, of which the tangent vector 
field is called left-invariant vector field on S . 

This fact motivates the idea that the left-invariant 
flows on the general Lie groups can be considered as gen- 
eral oscillators. In particular, the left-invariant flows on 
the non-Abelian Lie groups can be considered as non- 
Abelian oscillators. Let G be a Lie group and Y its Lie 
algebra. If v G T, then the left-invariant flow determined 
by v on the G is 



9) = 9 ■ ex P t v i * G R, g G G. 



(5) 



Different from those defined on the linear Euclidean 
spaces, the above dynamical system (5) is defined on the 
manifold G. 
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coordinate system on a manifold. Let {xx,x%, - ■ ■ ,x n ) 
be a local coordinate system about the identity e of Lie 
group G, n — dimG, and (0, • • • ,0) the coordinates of 
e. If (xi,x 2 , ■ ■ ■ ,x n ), (yi,V2, ■ ■■ ,Vn) and (zi,z 2 , ■ ■■ ,z n ) 
are the coordinates of <?i,<?2 £ G and 3132 £ G in this 
coordinate system respectively, then the multiplication 
function / = , / n ) can be defined as: 



21 = ■ ■■x n ;yi, ■■■Vn) 

%n fn (^1 1 %n I Vl 3 * 2/n ) ■ 



(G) 



Suppose that u € L, then the equations of the left- 
invariant flow determined by v are 



dt 



J2a4(x(t)), l<j< 



(7) 
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FIG. 1: (Color online) Order parameters mi (a) and mi (b) 
vs coupling strength K for different values of the rewiring 
probability p. All the simulation results are the average over 
100 independent runs corresponding to the case of n = 2. 



Next we consider the synchronization of network, of 
which each node is located a general oscillator on the 
Lie group G with its Lie algebra T. The collective syn- 
chronization behavior of the coupled general oscillators 
starts by randomly choosing an element v a ^ in F and 
an initial phase go € G for each node a. In the case of the 
general oscillator v a and go correspond to the frequency 
oj and c = e l9 in Eq. (3), respectively. 

It is worthwhile to note that, in general, the phase 
space of our general oscillator is both non-Abelian group 
in algebra, and non-linear space, which is a manifold, in 
geometry. Thus one may expect that there will be more 
interesting phenomena in this model. 

For a dynamical system on the manifold G, its equa- 
tions of motion should be written in a local coordinate 
system because, generally, there does not exist a global 



where (a^, a 2 , • • • a„) is the coordinates of v with respect 
to the basis • • • , ^-} x=0 , and I = (If ) is the 

matrix valued function: 



dyi 



y =o, l<i,j<n. 



(8) 



The set of equations (7) is the intrinsic dynamic of a 
single general oscillator. Note that, the multiplication 
functions are the local representation of the multiplica- 
tion in the Lie group, and the matrix valued function 
(8) is the map induced by multiplying by x at left. This 
map translates every vector in the Lie algebra to x. Eqs. 
(7) mean that the left-invariant flow determined by v are 
made from the vector field whose value at x is just the 
vector into which translating v. 

Similar to the case of coupled simple oscillators [l2^ . 
the coupled dynamics of general oscillators are: 



Ea?tf(* a (*))-£ E q 3 (x a (t) - x? (t)) 

i /3eA Q 
Xj , 



dx^jt) _ 
dt — 

a?(0) = 

(9) 

where qj is the coupling function assumed smooth, A a 
is the set of a's neighboring nodes, and x" is the co- 
ordinates of initial phase of node a. According to the 
existence and uniqueness theorem of ordinary differen- 
tial equations, Eq. (9) has an unique solution. Then, by 
the theory of Lie groups, this solution can be extended 
to an one-parameter subgroup, which corresponds to a 
left-invariant flow. 



III. SIMULATION RESULTS FOR SOME 
TYPICAL EXAMPLES 

In this section, we will show some analyses and sim- 
ulations on several typical examples. All the simula- 
tions are obtained based on the Watts-Strogatz (WS) 
network which can be constructed by starting with 
one-dimensional lattice and randomly moving one end- 
point of each edge with probability p. The network size 
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FIG. 2: (Color online) Order parameter m vs the coupling 
strength K. All the simulation results are obtained by aver- 
aging 100 independent runs. 



N = 1000 and average degree (k) = 6 are fixed. In 
our numerical simulation below, the coordinates a, are 
chosen randomly and uniformly in the range (—0.5,0.5), 
and the initial phases x" in (0, 2tt). All the numerical 
results are obtained by integrating the dynamical equa- 
tions using the Runge-Kutta method with step size 0.01. 
The order parameters (see below) are averaged over 2000 
time steps, excluding the former 2000 time steps, to allow 
for relaxation to a steady state. 



T n , the n-dimensional torus 



We firstly investigate the n-dimensional torus, T n 



S 1 xS 1 



x S , which is a connected compact Lie 



group. Denote G\ = T' n 
iR x iR x • • • x iR with L 1 



and its Lie algebra Fi = 
= (i.e. Abelian). This is 



a direct extension of the simple Kuramoto oscillators, 
and will degenerate to Karumoto oscillators when n = 1. 
It follows that an element in I\, that is an intrinsic fre- 
quency, is v — (iux, iu 2 , • • • , iu n ). 

It is easy to see that the general oscillation determined 
by v is periodic if u>i, u>2, • • • to n are integer-linear ly de- 
pendent, and quasi-periodic if wi, 0J2, • • *w n are integer- 
linearly independent. For the equations of motion we 
have 



d_ 

dt 



h 



/ sin(x" — Xi) 
\sin« -x&), 



(10) 



with the initial conditions x" (0) = x" , where K is 



the coupling strength. To characterize the synchronized 
states of the jth dimension, we use the order parameter 



1,2, 



,n, 



(11) 



where {•} signifies the time averaging. 

Fig. 1 displays the relations between the phase or- 
der parameter and the coupling strength for various val- 
ues of the rewiring probability p, where n — 2. Since 
the dynamical behaviors in different dimensions are in- 
dependent, and each one is equal to a simply Kuramoto 
oscillator, the synchronization behavior in each dimen- 
sion is the same as that of a Kuramoto oscillator [12j . In 
a word, the small-world effect ensures the global synchro- 
nization for sufficiently large coupling strength, and the 
higher disorder (i.e. larger p) will enhance the network 
synchronizability. 

This example can be used to describe the collective be- 
havior of the coupled double planar pendulums in a net- 
work. Furthermore, according to the Liouville- Arnold 
theorem, the phase space of a integrable system, if it 
is compact, is a n-torus. Therefore, this example can 
be widely applied to describe the synchronization phe- 
nomenon of integrable systems. 



B. GLo(3, R), the identity component of 
3-dimensional general linear group GL(3,R) 

GL (3,R), the identity component of GL(3,R), is a 
connected Lie group of dimension 9 (Actually, GL (3, R) 
is a open subset of R 9 ). Denote G2 = GLq(3, R), and T2 
its Lie algebra containing all the 3 x 3-real matrices. The 
coordinates (gij) for all g = (gij) 6 G2 can be addressed 
on the basis Eij of T2 (sec the case (a) of Appendix A), 
where 



1 
E n = I 





£33 = (000 
1 



In this circumstances the equations of motion are 

3 

l 7T a (t\a a . - 

h 

n-O: 



' lK,[l] E*?p(*Kj-ir £ ( 12 ) 



dt 



P =i 



/36A Q 



with the initial conditions xfj (0) = xfj , where (ag ) G T 2 
and K is the coupling strength. 

Since G2 is a subset of the Euclidean space R 9 , the 
coupling term of arbitrary pair of nodes a and (i can 



be directly written as xfj 
parameter is defined as: 



Therefore, the order 



N(N - 1) 



a<(3 



, (13) 
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The corresponding equations of motion are 




FIG. 3: (Color online) Order parameter m vs the coupling 
strength K. All the simulation results are obtained by aver- 
aging 100 independent runs. 



where x' 



(58) 



M 



and M 



max • 



{K(*)ll> 

In Fig. 2, we report the simulation results about the or- 
der parameter to as a function of K for different rewiring 
probabilities. Similar to the situation of T n , this system 
can approach to the global synchronized state for suffi- 
ciently large K when p is large enough. It is worthwhile 
to emphasize that the order parameter of different sys- 
tems can not compare to each other directly, since their 
dchnitions are not the same. 



C. SU(2), the special unitary group 

SU(2) = S 3 , the special unitary group, is a connected 
compact Lie group consisted of all anti-Hermitian trace- 
less matrices. Denote G3 = SU(2), and T3 its Lie algebra 
having the following basis 



Jl ~ \{ ) > J2 - \ ( -1 ) ' J3 - \ \ i 



Give an arbitrary element g <E SU(2), its local coordi- 
nates (xi, x 2 , X3) can be found from: 



Wi w 2 
-w 2 Wi 



and 



(expa;i J 1 )(expx 2 J2) (expa^Ja) , 

(14) 



W\W\ + W2W2 = 1. 



(15) 



dxf{t) 
dt 



K 

km 



E 

/3(EA Q 



sin -(*?(*) 



*?(*))> 



with the initial condition x"(0) 



'./0 



(16) 

j = 1, 2, 3, where 



if is the coupling strength. It is worthwhile to emphasize 
that the behaviors of this dynamical system is indepen- 
dent to the local coordinate systems. It follows that the 
local coordinates (14) in SU(2) can be replaced by 



g = exp (xiJi + X2J2 + X3J3) 



(17) 



which is convenient for the description of two-state quan- 
tum systems since the matrices Ji, J2, and J3 are relative 
to Pauli matrices as: 



-2iJ x = er 3 , 



-2iJ 2 



°2, 



~2iJ 3 



(18) 



Note that, a one-parameter subgroup U(t) C SU{2) 
can be considered as the time-evolution in the quantum 
mechanics. Since the subgroup e l< ^U(t), where ( is an ar- 
bitrary real number, represents the same time-evolution 
as U(t) does, we call two different nodes a and (3 are 
synchronous if they are only differ in a phase factor e 1 ^ , 
that is 











Therefore, if the equation 



a 13 



Wo W 



(eR. 
(19) 

(20) 



is hold for every 1 < a < (3 < N, the completely global 
synchronization is achieved. Consequently, for each pair 
of nodes a and (3, define the dynamical "distance" be- 
tween them as 



Huf3 



,fi 



(21) 



where | • | signifies the modulus of the complex number. 
Accordingly, the order parameter is 



r/>3 



= [(exp(-max{/i Q/3 |l < a < (3 < N}))] . 



(22) 



In Fig. 3, we report the simulation results about the or- 
der parameter to as a function of K for different rewiring 
probabilities. Analogously, this system can approach 
to the global synchronized state for sufficiently large K 
when p is large enough. Note that, e _1 < m < 1 for 
< I Ma/8 1 ^ 1- Thus in the uncoupling case (i.e. K = 0), 
the order parameter to approaches to 0.37. 

In terms of the Hopf bundle S 1 -> S 3 S 2 , the base 
space S 2 is the set of states of a two-state quantum system 
(see the ref. and Appendix B for details). Therefore, 
let a two-state quantum system be located at each node 
of a small-world network, the whole system can achieve 
the synchronized state if the coupling strength is large 
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enough. Note that, the special orthogonal group 50(3) 
has an isomorphic Lie algebra to T 3 for the existence of 
the double covering SU(2) — » 50(3), which is a homo- 
morphism. Therefore, the dynamical system 50(3) has 
the same synchronized behavior as that of SU(2). The 
case where each node is located with a top can be repre- 
sented by the system 50(3) [Tof. 

IV. CONCLUSION 

In this paper, we extended the concept of Kuramoto 
oscillator to the left-invariant flow on general Lie group, 
and studied the generalized phase synchronization on 
small- world networks. The left-invariant flow on Lie 
group can be used to represent many significant physical 
systems, such as the integrable systems (e.g. the double 
planar pendulums), the two-state quantum systems (e.g. 
Ising model), the motions of top, and so on. In particular, 
the dynamics of two-state quantum systems on complex 
networks are extensively studied recently 0,0,0, the 
present work can provide us a theoretical frame to inves- 
tigate their collective synchronized behaviors. 

It is intuitive that the low-dimensional systems are 
more easily to synchronize. However, one may note that 
in the regular networks (i.e. p = 0), the system T 2 
does not exhibit synchronized behavior while the ob- 
vious synchronization is observed in the higher dimen- 
sional systems GLq(3,R) and SU(2) when the coupling 
strength gets large. This is an interesting phenomenon 
that worths a further study in the future. 
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APPENDIX A: HOW TO CHOOSE THE 
COORDINATE SYSTEM 

In order to investigate the collective synchronization 
behaviors of the coupled non-Abelian oscillators, we have 
to integrate Eq. (9) numerically. Because of the using 
of local coordinate system, a question rises: Does the 
Eq. (9) express the behavior of the system in the whole 
manifold G? The answer is as follows. 

Case (a) Manifold G is a connected open subset of Eu- 
clidean space R n . In this case, we can directly use the 
natural coordinate system of R n as the local coordinate 
system of G. 

Case (b) G is a connected compact Lie group. In this 
case, we have exp(r) = G, so we can choose the coordi- 
nate system of T as the local coordinate system of G. 

APPENDIX B: RELEVANCE OF THE 
DYNAMICAL SYSTEM SU(2) 

According to the theory of fiber bundles, the condition 
(20) is equivalent to that g a (t) and g^it) is all in the 
same fiber of the Hopf bundle S 1 -> S 3 S 2 . In order 
to interpret the map p we consider the 2-sphere S 2 as the 
complex projective line CP 1 in which a point is written as 
[wi, W2], where w\ and W2 are two complex numbers and 
not all zero. Here, the point [101,102] is identified with 
[Xw\, XW2] for any nonzero complex number A. Then the 
map p can be written as 

p(wi,w 2 ) = [101,102]. (Bl) 

It is easy to show that the inverse image p^ 1 [101,102] is 
the sphere S 1 @- This fact tells us that, in the sense of 
our definition, the collective synchronization behaviors of 
the coupled non-Abelian oscillators SU (2) can represent 
the two-state quantum system. 



[1] S. H. Strogatz, SYNC-How the emerges from chaos in 
the universe, nature, and daily life (Hyperion, New York, 
2003). 

[2] Y. Kuramoto, Chemical Oscillations, Wave and Turbu- 
lence ( Springer- Verlag, Berlin, 1984). 

[3] A. Pikovsky, Synchronization (Cambridge University 
Press, Cambridge, 2001). 

[4] J. A. Acebron, L. L. Bonilla, C. J. P. Vicente, F. Ritort, 
and R. Spigler, Rev. Mod. Phys. 77, 137 (2005). 

[5] B. F. Schutz, Geometrical method of mathematical 
physics (Cambridge University Press, 1980). 

[6] T. Brocker, and T. T. Dieck, Representations of compact 
Lie groups (Springer- Verlag, 1985). 



[12] H. Hong, M. Y. Choi, and B. J. Kim, Phys. Rev. E 65, 
026139 (2002). 

[8] D. J. Watts, and S. H. Strogatz, Nature 393, 440 (1998). 
[9] T. Frankel, The Geometry of Physics, An Introduction 
(Cambridge University Press, Cambridge, 2004). 
[10] V. I. Arnold, Mathematical Methods of Classical Mechan- 
ics (Springer- Verlag, New York, 1978). 
[11] C. P. Herrero, Phys. Rev. E 65, 066110 (2002). 
[12] H. Hong, B. J. Kim, and M. Y. Choi, Phys. Rev. E 66, 
018101 (2002). 

[13] C. -P. Zhu, S. -J. Xiong, Y. -J. Tian, N. Li, and K. -S. 
Jiang, Phys. Rev. Lett. 92, 218702 (2004). 



